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Shape function in QED and bound muon decays 
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When a particle decays in an external field, the energy spectrum of the products is smeared. We 
derive an analytical expression for the shape function accounting for the motion of the decaying 
particle and the final state interactions. We apply our result to calculate the muonium decay 
spectrum and comment on applications to the muon bound in an atom. 


I. INTRODUCTION 

A bound particle decays differently than when it is 
free. Even in the ground state, due to the uncertainty 
principle, bound particles are in motion that causes a 
Doppler smearing of their decay products. In addition, 
if the charge responsible for the binding is conserved, 
daughter particles are subject to final state interactions. 

Binding effects partially cancel in the total decay width 
SHU- However, in some regions, the energy spectrum of 
the decay products can be significantly deformed. The 
range of the accessible energy can also be modified, by a 
participation of spectators. 

In this paper we focus on weakly bound systems in 
quantum electrodynamics (QED) where the bulk of the 
decay products remains in the energy range accessible 
also in the free decay. The slight but interesting redistri¬ 
bution in that region is governed by the so-called shape 
function S Here we present for the first time a 

simple analytical expression for S. 

The shape function was first introduced to describe 
heavy quarks decaying while bound by quantum chro¬ 
modynamics (QCD). It is employed in a factorization 
formalism based on the heavy-quark effective field theory 
(HQEFT) that separates the short-distance scale, related 
to the heavy-quark mass, from the long-distance nonper- 
turbative effects governed by the scale Aqq D , embodied 
in the shape function. In QCD it is a nonperturbative 
quantity that can be fitted using data but not yet derived 
theoretically. 

The shape function formalism has been defined also for 
quarkonium nu. Subsequently in [12] a quarkonium pro¬ 
duction shape function was obtained analytically. Ana¬ 
lytical results for the decay shape function in the’t Hooft 
model were obtained in Ca¬ 
in QED, the shape function has recently been com¬ 
puted numerically and applied to describe the decay of 
a muon bound in an atom m (so-called decay in orbit, 
DIO). The spectrum of decay electrons consists of the 
low-energy part up to about half the muon mass , and 
a (very suppressed) high-energy tail extending almost to 
the full m^. The shape function formalism applies only 
to the former, also known as the Michel region Ca¬ 
in this paper we will not be concerned with the high- 
energy tail. We note here only that it is also of great cur¬ 
rent interest because it will soon be precisely measured 
by COMET [16] and Mu2e [17] ■ The high-energy part of 
the DIO spectrum is a potentially dangerous background 


for the exotic muon-electron conversion search, the main 
goal of these experiments. That region has therefore re¬ 
cently been theoretically scrutinized cb na¬ 
il. FACTORIZATION IN MUONIUM 

The HQEFT is based on the heavy quark mass being 
much larger than the nonperturbative scale Aqcd- Sim¬ 
ilarly, in muonic bound states there exists a hierarchy 
of scales [20j: the mass of the decaying muon is much 
larger than the typical residual momentum, p. In 

a muonic atom we have 

p ~ m^Za, (1) 

while in muonium 

p ~ m e a, (2) 

where a « 1/137 is the fine structure constant and m e is 
the electron mass. 

With this observation, the factorization formula and 
the shape function for the muon DIO were derived in 
P3] using earlier QCD results [5j HJ [BjfTTIl . Here we fol¬ 
low an equivalent but a slightly more general approach 
fTJ to derive the differential rate for a heavy charged par¬ 
ticle decay in the presence of an external Coulomb field, 
neglecting radiative effects. We apply the result to find 
the decay spectrum of muonium. 

We concentrate on the muon decay p + —t e + u^v e but 
our results are general and apply to any QED bound state 
decay, provided that the momentum in the bound state 
is much smaller than the decaying particle mass. 

The decay amplitude is related to the imaginary part 
of the two-loop diagram depicted in Fig. [I] Integrating 



FIG. 1. Muon self-energy diagram whose imaginary part 
corresponds to the muon decay rate. Double line for charged 
particles indicates the electromagnetic interaction with the 
spectator electron that needs to be resummed. 
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over the relative neutrino momentum we express the dif¬ 
ferential decay rate as 

dr = 2Gfdm (h afj ) (3) 


where q is the sum of neutrino four-momenta and Gf is 
the Fermi constant |2T| [22]. The neutrino tensor is 


W pl/ = - 


K 2 
3(2tt)3 9 



q 2 ) ' 


(4) 


The charged-particle tensor h can be decomposed using 
five scalar functions that depend on q 2 and v-q = qo . Here 
v is the four-velocity of the bound state. In general, 

T h2Vf 1 Vi/ ih^e^ prj v p q 
+h±q p q 1 * + h 5 (q v v p + q p v u ) , (5) 

but since the neutrino tensor Q is symmetric under p -SA 
v, from now on we neglect the asymmetric part of h. 
Contracting the tensors and denoting Wi = Im(fy) we 
find that only two functions wi :2 suffice to describe the 
double differential spectrum, 


dr 

dq 2 dq 0 


G 2 f 

3(2tt) 4 



(? 2 - 9 o ) U, 2 




( 6 ) 

Functions Wi can be calculated in QED. Adopting 
Schwinger’s operator representation [23], we have instead 
of the free electron propagator 


1 

— m e 


f + f — m e 


(7) 


where is defined such that it does not contain any 
heavy degrees of freedom. The commutator of its com¬ 
ponents gives the electromagnetic field-strength tensor 
[7t m , tt v ] = —ieF p v where e is the muon charge. Formally, 


V = 2 (M 


Pin 


j + f — ■ 


(! - 75 )p 


M 


( 8 ) 


where | M) denotes the bound-muon state and k = m p v — 
q. Equation ^ is valid in the whole phase space. 

To simplify our considerations, we restrict ourselves to 
the Michel region where the electron is almost on-shell, 
k 2 ~ m p p and the time component of k is large, v-k p. 
This is the region where binding effects are most promi¬ 
nent. (Near the highest energies also the virtuality is 
much higher, k 2 ~ m 2 fll permitting a perturbative expan¬ 
sion of the electron propagator [H].) We neglect the elec¬ 
tron mass since the electron is highly relativistic [5J [23] . 
The only effect of the electron mass is an overall shift of 
the endpoint spectrum, just like in a free-muon decay. 

In the Michel region, the four-momentum k can be 
written as k = (v ■ k)n + <5fc, where n is a lightlike vector, 
n 2 = 0, and Sk ~ p. Neglecting terms suppressed by 


h pl , = 4 {2m li v p v v - v ■ kg pv - v„q p - v p q u ) 

1 


x ( M 


k 2 + 2 fy • n) (k ■ v) 


M 


(9) 


We cannot further expand the denominator since both 

T 2 

terms are of order m p p. We introduce A = — ; it will 

be useful to remember that A scales like the muon mo¬ 
mentum A ~ p ~ Zoc. We now define the shape function, 

5(A) = (M\6(X - n ■ n)\M) , (10) 

and obtain 


w i = 2nS(X) 
4 m 


2 m u 


w 2 = 


k ■ v 


^ttS(X) = M 


k ■ v 


w 1 . 


( 11 ) 

( 12 ) 


We have recovered the QCD scaling behavior [23]: func¬ 
tions Wi depend in the leading order only on the ratio of 
k 2 and v ■ k rather than on these two variables separately. 
Equation (10) reveals that the shape function is closely 


related to the momentum distribution of the muon in the 
bound state. However, due to gauge invariance we cannot 
replace n ■ n by the momentum in the ft direction. 


III. SHAPE FUNCTION 


Formula (10) is the same for muonium and for a muonic 


atom. Both systems are nonrelativistic, therefore the 
wave function needed to calculate the expectation value 
in (10) has the same analytical form. The only difference 


is its parameters and thus the physical scales that char¬ 
acterize the muon momentum p [see below, Eq. ©] . We 
now proceed to an explicit calculation of the function S 


in Eq. (101. 


The bound-state wave function follows from field the¬ 
ory via the Bethe-Salpeter equation [2B]. In the non¬ 
relativistic limit it reduces to the Schrodinger equation, 


P 

2/i 


V(r) i> s (r) = Eip s {r), 


(13) 


where p is the reduced mass of the system. In the case 
of a muonic atom, with the mass of the nucleus uin , 


m„mjv 

P = -:- « rn p . 

m p + m N 


(14) 


Subsequent formulas apply to muonium with the follow¬ 
ing substitutions, 


m p -A m e , 
m N -A m p , 

Z -A 1. 

For example, the reduced mass in muonium is 

m e nin 

p = ---« m e . 

m e + m p 


(15) 

(16) 


With this notation we also have p ~ Zap. 

As customary, Eq. (13) is written in the Coulomb 
gauge, with the electromagnetic four-potential given by 


eA p (x) = (V(r), 0,0,0), (17) 
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with V(r) = — — for a muonic atom or V{r) = —- for 
muonium. The determination of the shape function is 
especially convenient in the so-called light-cone gauge, 


n^A^x) = 0 . 


(18) 


In this gauge, the electron is effectively free up to effects 
quadratic in the electromagnetic field. The price for this 
simplification is a more complicated formula for the muon 
wave function. In the light-cone gauge, Eq. (101 takes a 
simple form in the momentum representation, 

S(A) = J (fc) S (A + n ■ tj ips (fc) • (19) 

We are neglecting terms of order ( Za ) 2 in the above ex¬ 
pression. To fulfill condition (18), we change the gauge, 


eA'Jx) = eA^x) + d^xix), 


with 


x(x) = X (x) = Za In (n • f+ r). 


( 21 ) 


This transformation changes the muon Schrodinger wave 
function in the IS state, ips(x), by an r-dependent phase 
factor, such that 

V’s( r ) V'(r) = ps(i") = (n ■ f + r)~ lZa 4>s( r )- 

( 22 ) 

After the transformation, the wave function is no longer 
rotationally invariant, since the gauge fixing distinguishes 
the direction of the outgoing electron. 

We use the Schwinger parametrization, 


r(q) 

A a 


POO 

/ dtt a ~ l exp [— At] 

Jo 


(23) 


(20) to Fourier-transform Eq. (22) 


ip{k) = J d 3 r exp (-ife-f)^(r) 


iZa 


8^/ p 3 Z 3 a 3/ K 3 ( p 2 Z 2 a 2 + k 2 


iZa. 


r(iZa) sin (inZa) ^ 2 % 2 a 2 _|_ £ 2 ^ y2 iZa — in ■ itj J 2 (^xZa — in ■ k^j 


p 2 Z 2 a 2 + k 2 


— fi(i + Za) 


(24) 


We integrate in (19) first over n ■ k using the delta-function, then over k±, components of k perpendicular to n, 


_ 2/J 3 Z 6 a 6 3A" + 6 A/1 + /r 2 (4 + Z 2 a 2 ) arctanf 

“ 3sinh(7rZa) [A 2 + p 2 Z 2 a 2 } 3 


(25) 


The exponential function in (|25|) arises from 
(pZa + iX) 


Further, Eq. (25) allows us to better control the expan- 


— iZa 


appearing after integrating \xp(k)\ 2 

with the delta-function in The leading behavior 

can be understood with the help of integral 

1 1 


d 2 k _l 


(k 2 ± + A 2 + p 2 Z 2 a 2 


(A 2 + p 2 Z 2 a 2 ) 3 ' 


sion and the resummation of the p effects in the decay 
spectrum. This cannot be done so easily with a numeri¬ 
cal evaluation |14j , as is especially clear when we analyze 
the first three moments, useful in HQEFT for constrain¬ 
ing possible forms of the shape function. 

The zeroth order moment gives just the normalization. 
With the normalized wave function in Eq. (19), the shape 


(26) functions (251 and (27) are automatically normalized to 


The result (25) contains subleading terms, related to the unity; this is a consequence of the definition (101. 


Coulomb potential in (101, required by the gauge invari¬ 


ance. At the current stage of calculations 5(A) is explic¬ 
itly gauge independent. We can drop subleading terms 
to obtain a leading-order formula, 


When the subleading terms are neglected the first mo¬ 
ment of the shape function (27) vanishes, 


5(A) = 


8 p 5 Z 5 


J dXXS(X) = {n- tt) = 0 + 0(Z 2 a 2 ). 


(28) 


cr 


(27) 


37t [A 2 + p 2 Z 2 a 2 ] 

The analytical formula obtained here is useful for sev¬ 
eral reasons. First of all, counting powers and remem¬ 
bering that A ~ p, we find that 5(A) ~ - ~ This 
reminds us that 5(A) is a nonperturbative object and 
explains why its effect on the spectrum can be quite dra¬ 


matic, as we shall see in muonium in Sec. IV 


Naive power counting in the left-hand side suggests a re¬ 
sult linear in Za. That leading part vanishes, similarly 
to the first moment of the B-meson shape function. A 
contribution linear in -2- ~ Za is absent due to the 
CGG/BUV theorem [3JII]. Moments of the shape func¬ 
tion are related to matrix elements of local operators in 
the heavy particle effective theory. Operators of dimen¬ 
sions 3 and 5 exist. A dimension 4 operator that could 
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generate, in the leading order, a nonvanishing first mo¬ 
ment is missing. A nonzero first moment can only appear 
at the subleading order mm- 

The second moment is related to the square of the 
average momentum in the direction of ft, 

I dXX 2 S(X) = | (p 2 ) + O (Z 4 a 4 ) . (29) 


In contrast to the first moment, there is no cancellation 
here and the naive counting correctly predicts a result 
quadratic in Za. Therefore we do not need subleading 
corrections in Eq. © to calculate ( [29| . This moment 
characterizes the width <j\ of the region smeared due to 
the shape function effects. As expected, it is of the same 
order as p: a\ = This is similar to the HQEFT 

where the second moment is also related to the average 
kinetic energy of the heavy quark inside a meson. 

In muonic aluminum, the stopping target of the 
planned conversion searches (Mu2e and COMET), the 
shape function effect is sizeable since o\ ~ 6MeV, and 
has been precisely measured by TWIST m ■ In the case 
of the muonium the effect is much smaller, a\ ~ 2 keV, 
and is negligible except near the end of the spectrum. 

In Fig.[2]we plot the shape function for Za = 0.25. The 
width is proportional to p , suggesting that the dominant 
effect is due to the muon motion in the initial state. 



FIG. 2. Leading shape (271 function for fj, = 1 MeV and 
Za = 0.25. The width is a\ = /u^|. To illustrate it better 
we included two vertical dashed lines at the ends of ax region. 


Finally, we would like to point out that the formula 


(271 can guide phenomenological models of the shape 
function in QCD. Some QCD bound states can be de¬ 
scribed with the help of effective theories similar to non- 
relativistic QED [28j[29]. For example, Ref. [30] postu¬ 
lated a similar functional form of the shape function, 


S( A) = N {X X l 1 b) 2 + a 2 e W 0 O- - A )> ( 3 °) 


with parameters a, b to be fitted from data. Our function 
has a higher power of the denominator, therefore does 
not require an artificial restriction of its support by 6 
functions, because its tails are sufficiently suppressed. 


IV. MUONIUM SPECTRUM 

Having obtained the shape function, we can calcu¬ 
late the muonium spectrum using ([6]). After an inte¬ 
gration over q 2 , the electron energy is given by E e = 
— <7o + 0(Z 2 a 2 ). The shape function formalism can 
be interpreted as a replacement of the zero-width on- 
shell relation for the electron by a finite-width shape 
function 5(A). (If 5(A) in the functions Wi is replaced 
by the Dirac-delta on-shell condition, the free-muon de¬ 
cay spectrum results.) Since the muon is almost at rest, 
the smearing is negligible far from the free muon decay 
endpoint, the only region where the spectrum is quickly 
varying with the electron energy. 

We ignore the tail of the spectrum at energies higher 
than the free endpoint plus several am e . It is very sup¬ 
pressed and its evaluation requires perturbative correc¬ 
tions due to hard photons m ns. We also ignore the 
lowest region of the spectrum where positronium can be 
formed m- 

For illustration, Fig.[3]shows the muonium decay spec¬ 
trum in the vicinity of the free muon endpoint. The 
extent of the region affected by the shape function corre¬ 
sponds to the smearing width cr \, denoted by two vertical 
lines. In this region the slope of the spectrum is propor¬ 
tional to the shape function 5(A) and therefore is of the 
order of - ~ —. 

P <T\ 



FIG. 3. Endpoint region of the muonium decay. Electron en¬ 
ergy is parametrized in terms of e, such that e = 10 4 Ee ~ fimax , 

“max 

m 2 +m 2 

E max = —-. Dashed (solid) line shows the free-muon 

(muonium) spectrum. Vertical dotted lines emphasize the 
size of the region that is smeared due to binding effects. 
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Note that the free-muon decay, denoted with the 
dashed line, resembles a step function. This is an arte¬ 
fact of the very narrow width of the region shown in this 
figure. In fact, the free-decay spectrum varies with e to 
the left of the step and vanishes to the right of it. 

V. CONCLUSIONS 

We have derived an analytical formula for the shape 
function and used it to calculate the muonium decay 
spectrum. Shape function moments were analyzed and 


compared with appropriate expressions in HQEFT. Our 
analytical formula may also have a limited application to 
describe nonrelativistic QCD systems. For now, the ana¬ 
lytical expression for the shape function has improved our 
understanding of the approximations used in the deriva¬ 
tion of the muon DIO spectrum. 
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